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ABSTRACT 



o 



The study of ionized gas morphology and kinematics in nine eXtremely Metal- 
Deficient (XMD) galaxies with the scanning Fabry-Perot interferometer on the SAO 
_ V , 6-m telescope is presented. Some of these very rare objects (with currently known 

■ range of 0/H of 7.12 < 12+log(0/H) <7.65, or Zq/35 < Z <Zq/10) are believed to 
(~| \ be the best proxies of 'young' low-mass galaxies in the high-redshift Universe. One 
pi |. of the main goals of this study is to look for possible evidence of star formation 

(SF) activity induced by external perturbations. Recent results from HI mapping of 
a small subsample of XMD star-forming galaxies provided confident evidence for the 
^ important role of interaction-induced SF. Our observations provide complementary or 

' new information that the great majority of the studied XMD dwarfs have strongly 

disturbed gas morphology and kinematics or the presence of detached components. 
We approximate the observed velocity fields by simple models of a rotating tilted thin 
disc, which allow us the robust detection of non-circular gas motions. These data, in 
turn, indicate the important role of current/recent interactions and mergers in the 
observed enhanced star formation. As a by-product of our observations, we obtained 
\ data for two LSB dwarf galaxies: Anon J012544-I-075957 that is a companion of the 

■ merger system UGC 993, and SAO 0822-1-3545 which shows off-centre, asymmetric, 
• ' low SFR star-forming regions, likely induced by the interaction with the companion 

g ; XMD dwarf HS 0822-^3542. 

. Key words: galaxies: dwarf - galaxies: interactions - galaxies: kinematics - galaxies: 

evolution - galaxies: star formation - galaxies: individual: (UGC 772, HS 0122-1-0743, 
SBS 0335-052E and W, HS 0822-^3542, SDSS J1044-F0353, SBS 1116+517, 
^ ^ . SBS 1159-1-545, HS 2236-1-1344, SAO 0822+3545) 

X 

1 INTRODUCTION to intensive searches for new XMD galaxies in sev- 



J r 1 • • 1 era! recent large surveys, fe.g.. iMelbourne fc Salzeij I2OO2I: 

the group ot very rare dwart galaxies with gas metal- [77 ^ i-I-, A, ., — 7 — r"r,fi^o ' 1 — T 

,. . . ^ _ , , „ i.,r. , , ,TTN „ ^.r- lUgrvumov ot al. 2003; Pustilnik ot al. 2003c: Kniazov et al 

licities Z below Zp/lO ( o r 12+ log(0/H <7.65, e.g., re- „ ., , , „„„r „ 1 p 11 onn^ 

177 \ — p7t — r~r JTT^T^N , \r , i2003; Pustilnik ct al. 2005a: Brown, Kcwloy & GoUcr 200< 

view by i Kunth fc Osthn| ((20001)), are called eXtremely r^Tl (^nn4 lr:,,.„,„ .1 1 b^^^^ „rv,„r,. ^tin«; 



4- va^t^ i • q t, u- 4- ■ IPapaderos et al.1 120081: iGuseva et all |2009|. among others, 

Metal-Dehcient, or AMD galaxies, such obiects comprise , „ ■ \ i r j_i i i 

° and references therein), the number ot currently known such 

only about two per cent ol known blue compact galax- i • , . i i tt • • i . mi ■ i 

. objects m the local Universe is about 100. there is also 

les (BCGs), low-mass galaxies with active star formation . . , ,. . i i ■ ^ j_ i i i -n^ 

^ I "-r 1 1 — " — I progr ess m identihcation ot such objects at large redshilts 

(e g. [ Pustilnik et al. | | 2005a| ) . The latter in turn, comprise tKakazu. Cowie fc HuB . 
about 6 per cent ot th e entire local dwarl galaxy popu- 
lation l|Lee et al.l |2009| , and references therein). Thanks The interest in these atypical, local Universe galaxies 

is caused by the similarity of their properties to low- 
mass galaxies in the early Universe, when the typical 
* Based on observations obtained with the Special Astrophysical metallicity of galaxies was one-two orders of magnitude 

Observatory RAS 6-m telescope. smaller than the solar one. Therefore, understanding the 

t moisav(agmail.com(AVM),sap@sao.ru{SAP),akniazev@saao.ac.za(AYS{)ecific processes of star formation (SF) and evolution 
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in nearby galaxies with such low metallicity will pro- 
vide insights into similar processes in young galaxies at 
high redshifts. To better understand various aspects of 
XMD galaxies star formation and evolution, we conduct 
an extensive study of a subsample of X MD galaxies in 
optica l, NIR and radio domains fe.g.. | Pustilnik at al 
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One of the fundamental questions of SF in low-mass, 
gas-rich galaxies is the trigger mechanisms of starbursts: 
whether the intrinsic instabilities in gas discs are capa- 
ble to produce the observed level of enhanced (bursting ) 
SF (e.g., models by Pelupessy . van der Werf. fc Ick3 (|2004i ) 
and IPi Matteo et al.l I 2OO8); ob serva t ions of Virgo cluster 
dirs by Brosch. Almoznino. fc Ifelleij (|2004) and references 
therein) , or external mechanisms (such as tidal disturbances 
and mergers) are necessary in a significant fraction of ob- 
served starbursts. 

From the analysis of the density distribution in 
BCGs in compa rison to more typical late-type dwarfs, 
ISalzer fc NortonI Cl999) concluded that the former have a 
higher gas concentration and, therefore, BCGs represent a 
special group of galaxies, more susceptible to internal insta- 
bilities and related starbursts. But, it is unknown whether 
this is an inherent property of BCG progenitors, or this 
higher gas concentration appeared due to a minor merger, 
or as a result of a recent interaction with a barely detectable 
companion. If due to mergers or interactions, BCG progen- 
itors could be the quite common late-type dwarfs. To dis- 
entangle the possible scenarios one needs a large, represen- 
tative sample of starbursting galaxies with well studied and 
understood kinematics of both HI and ionized gas. This is 
required to determine the prevalence of external or internal 
disturbances. 

First statistical indications for the possible importance 
of interaction s to trigger starburst s in BCG progenitors was 
presented by iTavlor et al.l (|l995l ) in their search for HI- 
rich companions in a small BCG sample. For a larger BCG 
sample, the conclusion on the h igh fraction of i nteract ion- 
indu ced starbursts was m ade bv lPustilnik et akl l|2001al . see 
also iNoeske et all (|200ll )). based on the comparison of ob- 
served relative distances and 'threshold tidal distances', at 
which a perturber galaxy should induce shocks in the gas 
disc and pro duce dissip ation of gas angular momentum (as 
suggested bv llckdll985l ). While that study was indicative of 
the important role of collisions in starbursts of gas-rich low- 
mass galaxies, more elaborate tests should be used to clear 
up the real role of external factors in starbursts. In turn, 
numerous studies of galaxy morphology and k i nematics at 
redsh ifts of up to ~1 (e.g., iPuech et al] I2OO6I : lYang et ahl 
|2008| ) show clear indications on the importance of external 
mechanisms on the enhanced SF in a significantly larger 
galaxy fraction than at the current epoch. 

In many cases, the traces of interactions can be hidden 
if only images of stellar components are available. There- 
fore, gas morphology and kinematics provide a better tracer 
of interaction-induced starbursts. To this end, HI mapping 
is a good way for galaxies with sufficiently large HI flux. For 
less gas-rich and/or more distant galaxies, the study of the 



ionized gas kinematics can be a good, complementary tech- 
nique. Besides, this method provides usually a finer angular 
resolution and allows comparison of the global gas kinemat- 
ics derived from HI mapping with the gas motions within 
the region of higher density, around the s i tes of SF. T hese 
kind of studies, initiated by lOstlin et al] l| 19991 . 1200 ll ) for 
a small sample of luminous BCGs, already provide good 
evidence for merger-induced SF in such objects. Some in- 
dications for the sam e phen omenon have been observed by 
iGarcfa-Lorenzo et all ||2008|) for a group of five other lu- 
minous BCGs and by|P erez- GaUego et all (|2009l ) for repre- 
sentatives of the similar group of Luminous Compact Blue 
Galaxies. 

The starbursting XMD galaxies are especially inter- 
esting in this aspect, since the results may have implica- 
tions for cosmological young galaxies in the high-redshift 
Universe. The evidence of significantly disturbed optical 
morphology and interactions, indicating an unrelaxed state 
of many XMD galaxies, was already noticed by us (e.g., 
'Pustilnik ct al.' 2004b, '2 005bl . l2006 h and more recently by 
,Papadcros ct al. (2008). How does this relate to the origin 
of XMD galaxies? 

There are several options for evolutionary scenarios of 
XMD dwarf galaxies, alread y emphasized in the literature 
(e.g., IPustilnik etlU l2005bh . Two of them relate to gas- 
rich discs, which are very stable locally, and being isolated, 
have been evolving very slowly (like LSB ga laxies). Some o f 
them might even remain dark objects (e.g., |Pustilnil3l2008l ) 
(and, thus, 'pristine' like protogalaxies at epochs of the ma- 
jor galaxy formation). Such galaxies would be very difficult 
to discover/recognise before they (by chance) have experi- 
enced a strong external perturbation due to an interaction 
with a galaxy-sized mass. When this happens, such galaxies 
can lose their global stability and experience a substantial 
SF burst. Due to their slow evolution, during the starburst 
they would appear as gas-rich and very metal-poor (XMD) 
objects. Therefore, the study of gas kinematics for the sam- 
ple of XMD galaxies can shed light on their origin as a group, 
or can give clues to the diversity of their properties and their 
finer classification. 

In this paper we present the results of an Ha study of 
the ionized gas kinematics in nine XMD starbursting galax- 
ies (BCGs) conducted with the SAO 6-m telescope's scan- 
ning Fabry-Perot Interferometer (FPI). These galaxies con- 
stitute a part of a larger XMD galaxy sample, for which we 
conduct a multiwavelength study of their properties. Most 
were selected for this FPI Hq study due to their unusual or 
disturbed optical morphology. For a fraction of them, the 
HI mapping data are available, which provide an opportu- 
nity of the combined view on the gas kinematics. We plan to 
continue this project for a larger sample of XMD starburst- 
ing galaxies in order to improve statistics and to search for 
possible differences in the properties of the ionized gas kine- 
matics. 

Due to volume limitations, the scope of this paper is 
kept mainly to observational results for the studied XMD 
galaxy subsample and their preliminary analysis. The more 
advanced analysis, involving multiwavele ngth data and new 
methods of modelling (like suggested bv lBarnes fc HibbardI 
200!|), will be presented for individual galaxies in forthcom- 
ing papers. This paper is organised as follows: in Sect. [2] we 
describe observations and their reduction. Sect. [3] presents 
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Table 1. Log of the 6-m telescope FPI observations 



Object Date Exposure ang. rcsol. 







time (s) 


(arcsec) 


(1) 


(2) 


(3) 


(4) 


SDSS J0113+0052 


11.01.07 


36x200 


4.0 


HS 0122+0743 


10.01.07 


2x36x90 


2.3 


SBS 0335-052W 


14.01.08 


2x36x150 


1.4 


SBS 0335-052E 


14.01.08 


2x36x150 


1.4 


HS 0822+3542 


10.01.07 


36x250 


1.4 


SDSS J1044+0353 


14.01.08 


36x150 


1.4 


SBS 1116+517 


10.01.07 


36x250 


1.5 


SBS 1159+545 


14.01.08 


36x170 


1.8 


HS 2236+1344 


08.09.05 


36x250 


1.0 



well (see discussion in lMoiseev fc Egorovll200^ '). The profile 
fitting results were used to construct the two-dimensional 
field of line-of-sight velocities of ionized gas, the map of ve- 
locities dispersion, and images of the galaxy in the Ha line 
and in the neighbouring continuum. We estimate the ac- 
curacy of line-of-sight velocity measurements through the 
m easured S/N ratios , using the relations given in Fig. 5 
of Im oiseev fc Egorovl l|2008l ). The accuracy was found to 
be 1-5 km in the bright HII regions and amounts 
to 10-15 km s~^ in the regions with the minimal surface 
brightness, where the emission line shows up at the S/N 
ratio level of S/N « 5. In Fig. [T] we show the examples 
of representative Ha-line profiles for three program galaxies 
SBS 0335-052W and E and SBS 1116+517. 



the description of data analysis and the summary of results. 
In Sect.|4]the kinematic properties are discussed in more de- 
tail and are compared with other data on these BCGs, and 
their likely interpretation is suggested. Sect. [5] summarises 
conclusions of this study. 



2 OBSERVATIONS AND DATA REDUCTION 

The observations of XMD galaxies were conducted with 
the m ultimode instrument SCORPIO ( Afanasiov & Moiseev 
|2005[ ) installed in the prime focus of the SAO 6-m telescope 
(BTA) during 3 runs (see Table [1] for details and Table [2] for 
galaxy main parameters). The desired spectral interval in 
the neighbourhood of the redshifted Ha line was cut using a 
set of narrow-band filters {FWHM = 15 - 20A). The width 
of the free spectral interval between the neighbouring orders 
of interference was equal to 13 A (about 600 km s~^). The 
resolution of the interferometer {FWHM of the instrumen- 
tal profile) was 0.8 A (35 km s'^) for a 0.36A/channel sam- 
pling. The detector was a 2048x2048 EEV 42-40 CCD oper- 
ating in the 2x2 and 4x4 (for SDSS J0113+0052=UGC 772 
and HS 0122+0743) pixel instrumental averaging mode in 
order to reduce the readout time. The corresponding im- 
age sampling was 0'.'36 and 0'.'72 pixel" ^ respectively. The 
resulting field of view has a size of 6'.1. In column 4 we 
also present the final angular resolution of our data which 
is a combination of the real seeing and the data reduction 
smoothing. 

We took a total of 36 successive interferograms per ob- 
ject with different gaps between the plates of the FPI. The 
observational data were reduced using an IDL-based soft- 
ware package l|Moiseevll2o"o3 : iMoiseev fc Egorovl l2008f) . Af- 
ter primary reduction, night-sky line subtraction, and the 
wavelength calibration, the observations were reduced to the 
form of a data cube, where each pixel in the field of view 
contains a 36-channel spectrum. We removed the ghost im- 
ages that appear in the plat es of the inte rferometer using 
the algorithms described bv iMoiseev &: Eg orov ( 2o3^ ). To 
check the quality of ghosts removal on the data of extended 
objects, we observed HS 0122+0743 and SBS 0335-052E,W 
successively in two fields turned in the position angle. 

We fitted the Ha emission profiles to the Voigt function, 
which in most cases describes the observed contour fairly 



3 RESULTS AND ANALYSIS 

In Table [2] we present the summary of the main observa- 
tional parameters of the studied XMD galaxies. Since the 
accuracy of our total Ha flux estimates (10-30 per cent) 
is in general worse than that for dedicated Hq photome- 
try, we give our values only for unknown Ha fluxes. For 
the remaining objects, we cite data from the literature. Ab- 
solute magnitudes, Mb, are calculated for distances taken 
from NElJI, except for HS 0822+354 2 and its c ompanion 
LSB dwarf. As discussed in iPustilnik et al.l (|2010i ). they are 
situated in the region with large peculiar negative radial ve- 
locitie s (~300 km s~^) due to the effect of the huge Local 
Void (|Tu11v et al.l bOOSl ). So to derive their distances, this 
peculiar velocity was accounted for. 

The dominant or important mode of gas motions in late- 
type, low-mass isolated galax;ies (even for objects with lumi- 
nosities corresponding to My fainter than -14.0) is rotation 
in a gravitational field of a Dark Matter (DM) halo and 
baryonic disc (e.g.. lBegum et aL I l2006l .[ 20081 . and references 
therein). Two types of processes can significantly affect the 
regular velocity field attributed to rotation: the star forma- 
tion activity and tidal interactions. 



3.1 Kinematics of SF induced shells 

The strong SF activity releases a significant amount of ki- 
netic and thermal energy in short periods, and it results 
in the formation of hot bubbles and shells with typical 
sizes of a hundred pc to about one kpc. Such shells, with 
maximal observed line-of-sight velocities in the range of 
ten to a hundred km s~^ are recognised in both long-slit 
spect ra and 2D Ha velocity fields of star-forming galaxies 
(e.g.. lMartinlll996l. Il998l: IZasov et al.ll2000l: [Pustilnik et al l 



l2003b':'Lozinskava et alll2006l : lMartinez-Delgado et al.ll2o'o7[ 

Bordalo ct al. 2009, among others). In principle, the large- 
scale velocity fields in such expanding ionized shells are suf- 
ficiently simple. However, in case of dwarf galaxies, which 
have relatively small rotation velocities, shells can signifi- 
cantly perturb the disc velocity field. In Appendix [Xj we 
present the simulated velocity field of dwarf galaxies with 
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shells for several simple cases. Using the characteristic pa- 
rameters typical of our sample (disc size, inclination, rota- 
tion curve, etc.), we describe the technique for such shells 
recognition and masking in the observed velocity fields. For- 
tunately, in many simple theoretical scenarios it is possible 
to extract the real parameters of global rotation, even when 
the shell expansion velocity is larger than the line-of-sight 
projection of global rotation on given radii. 



3.2 External perturbations of velocity fields 

It is not yet well understood how the strength of starbursts 
and that of related expanding bubbles are connected with 
internal or external triggers. But, since the external pertur- 
bations (e.g., due to galaxy collisions or extragalactic cloud 
infall) provide disturbances which are additional to internal 
ones, they are naturally expected to trigger the stronger dis- 
turbances in galactic gas and to result in more intense star 
formation. Besides, in general, the characteristic timescales 
of singular star-forming events in isolated low-mass galax- 
ies (of ten to several tens Myr) are shorter than those for 
the events induced by external interactions. The latter are 
comparable to rotation periods, and usually are of the or- 
der of a hundred Myr and more, the typical timescale for 
galaxy response to a collision- induced perturbation (e.g., 
iDi Matteo et al.]|2008l . and references therein). 

This is due to different spatial scales involved in such SF 
episodes. For internal perturbations we usually have some 
local enhancement, while for external ones the gas flows 
on scales of the whole galactic disc are involved. Therefore, 
one can expect rather different effects in isolated and in- 
teracting starbursting galaxies. In the former, the gas ve- 
locity field should display the characteristic rotation pat- 
tern with rather localised disturbances due to the starburst- 
related shells. More or less typical dirr galaxies can serve 
as templates for such kinematics. In the interaction-induced 
starbursts one can observe, besides the overall rotation and 
the appearance of shell kinematics, the effects of tidally dis- 
turbed velocity fields to various degrees. In extreme cases of 
a major merger, the strong disturbance of the initial equi- 
librium rotation field can almost wash out the regular com- 
ponent. 

Apart from non-destructing (fly-by) collisions, various 
types of mergers (major, intermediate mass-ratio, and mi- 
nor [Mi:M2 > 5]) can trigger late-type galaxy starbursts. 
Depending on the stage and mass proportion, the veloc- 
ity fields can show various patterns - from very disturbed 
with no traces of regular rotation to more or less regular 
rotation with residual motions (e.g. counter-rotation in the 
central region, where a smaller component has sunk as the 
result of minor merger). As recent simulations d emonstrate 
fe.g.. lKronberger et allbooel :! Pedrosa et al bops '), the bifur- 
cations in velocity curves should appear (that is prominent 
asymmetries in receding and approaching branches of rota- 
tion curves) as a result of close recent encounters. Various 
appearances of complex velocity fields resultin g from major 
merge rs of disc galaxies were demonstrated bv lJesseit et al.l 



3.3 Fitting results 

From the above discourse it is clear that it is quite natu- 
ral to analyse Ha velocity fields in the studied XMD dwarf 
galaxies by trying to fit them with a model of a fiat rotating 
disc and to compare the residual velocity field against vari- 
ous patterns characteristic of shells, warps, central counter- 
rotation, discontinuities, etc. Of course, the morphology of 
both stellar continuum and Ha emission should be taken 
into account in the course of the results' interpretation. 

We analysed the velocity fields using the well-known ap- 
proximation of a thin rotating disc ('tilted-rings' method), 
which is generally accepted in analysis of gas kinematics. 
For the method description an d references to the original 
works, see [m oiseev et al.l(|2004l ). This method calculates the 
disc orientation parameters (major axis position angle PA, 
inclination i, systemic velocity Vsya) together with radial 
variations of the rotation velocity Vrot- In special cases of 
non-circular motions, we can also estimate radial dependen- 
cies of PA and Vays- 

In Table |3]the 'best-fitting' model parameters are sum- 
marised: systemic velocity Vsys, position angle PA and incli- 
nation angle i, along with relevant brief comments. Unfortu- 
nately, from the kinematic model we could only derive i for 
three galaxies, whose inclinations are indicated with errors 
in the Table O For J0113-f 0052 and SAO 0822-1-3545, we 
adopted the inclinations from published HI maps. In the re- 
maining cases, to estimate the inclination angle we adopted 
the approach based on the morphology of the overall Ha 
emission and the observed ratio of minor to major axis of 
Ha 'ellipsoid' p=a/b. Assuming that this refiects the real gas 
distribution for an inclined 'thick' (typical of dwarf galaxies) 
disc with the intrinsic axial ratio q, we estimated the inclina- 
tion i with well-known formula [cos{i)]'^ = {p^ — q^)/{l — q^). 
The estimate of parameter q was adopted according t o 
formula A6 from lStavelev-Smith. Davies fc KinmanI (| 19921 ). 
We understand that in many cases, where we infer strong 
disturbances in morphology and velocity field, the above as- 
sumption can be hardly justified, and hence the adopted 
inclination angle and the related amplitude of rotation ve- 
locity can be rather uncertain. But the concrete choice of i 
does not affect the overall radial dependence of the model ro- 
tation velocity. A more detailed description of the data and 
their model fits, along with discussion of relevant kinematics 
and morphology in Hl-line and other points are presented 
for each individual galaxy in Section |4l 

For each of the nine program XMD galaxies and for 
two 'by-product' LSB dwarf galaxies, we show the related 
graphic materials in Figures I2I4I The colour images of ve- 
locity field data are available in electronic version of the 
journal. They are arranged in columns of 6 images in the 
following order. The first (top) panel displays continuum 
image in ^-filter from the SDSS or in the blue band from 
the digitised Palomar Observatory Sky Survey (POSS). 
For SBS 0335-052E W, the BTA B-band images are used 
from IPustilnik etall l|2004bl ') . The brightness distribution is 
shown in logarithmic scale. The second panel shows the im- 
age of integrated emission in the Ha line in logarithmic scale. 
The third panel shows the velocity field, both in colour and 
by isovelocity lines. In the forth panel, we show the 'best- 
fitting' model of thin tilted rotating disc. The fifth panel 
presents, in colour, the residual velocity field (observed mi- 
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Table 2. Parameters of nine XMD and two by-product (below the line) galaxies 



lAU name Coord. (2000.0) Vhdl BtotJ M^* O/Ht F(Ha)§ D§§ Scale Alternative 

R.A. Dec. km s^^ mag mag Mpe pcarcsec"^ name 

h m s o I II 



(1) (2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(11) 


SDSS J0113+0052 01 13 39.40 


+00 52 27.9 


1176+ 1 


16.28^ 


-14.70 


7.24 


6.8 


16.3 


79 


UGC 772 


HS 0122+0743 01 25 34.18 


+07 59 22.2 


2940+ 2 


15. 48^ 


-17.76 


7.63 


20.0^ 


40.3 


195 


UGC 993 


SBS 0335-052W 03 37 38.46 


-05 02 36.3 


4038+ 2 


19. U'' 


-14.73 


7.12 


1.44 


53.8 


261 




SBS 0335-052E 03 37 44.04 


-05 02 37.6 


4053+ 2 


16.95'! 


-16.92 


7.29 


32.0-* 


53.8 


261 




HS 0822+3542 08 25 55.43 


+35 32 31.9 


727+ 2 


17.92^ 


-12.49 


7.45 


7.6'^ 


13.5 


65 




SDSS J1044+0353 10 44 57.84 


+03 53 13.2 


3865+ 1 


17.625 


-16.20 


7.44 


14.5 


53.8 


261 




SBS 1116+517 11 19 34.29 


+51 30 10.7 


1342+ 2 


17. 14^ 


-14.70 


7.51 


15.2 


23.1 


112 


Arp's 


SBS 1159+545 12 02 02.36 


+54 15 50.1 


3592+ 1 


19. 00^ 


-14.60 


7.49 


5.5 


52.2 


255 




HS 2236+1344 22 38 31.15 


+14 00 28.6 


6162+ 2 


17.88^ 


-16.90 


7.50 


12.0^ 


86.4 


419 




Anon J012544+075957 01 25 44.18 
SAO 0822+3545 08 26 05.59 


+07 59 57.0 
+35 35 25.7 


3012+ 3 
742+ 2 


17.0:2 
17.56^ 


-16.30 
-12.85 




1.5 
0.4 


40.3 
13.5 


195 
65 





(II) " systemic velocity on our Ha data. For UGC 993 this is a mean of Vsys for E and W components; 

for SAO 0822+3545 - this is HI velocity from Chengalur et al. (2006) 
{X) — Photometry data are from: ^ ^mokerj_^avies_^^xon (1996); ^ Pustilnik et al. in prep.; ^ Pustilnik et al. (2003a); 

^ Pustilnik et al. (2004b); ^ NED/SDSS transformed from g,r to B, ^ Moustakas fc Kennicutt (2006*) for Ha flux, 
(t) - in units 12+log(0/H). Data are from Ugrvumov et al. (2003) . Kniazev et al. (2003) . Izotov fc Thuan (2007) . 

Izotov. Thuan fc Guseva (2005) . 
(*) — corrected for A^ (according to ^cUegelj_^inkbeiner_^^ouglas (1998)), with distances from Col. 9 
(§) - in units of lO"^** erg cm"^ Ha fluxes adopted from literature are coded as photometry data. 

(§§) - from NED [Virgo-inflall; Ho=73], except HS 0822+3542/SAO 0822+3545, for which additional correction is adopted (see text). 



Table 3. Parameters of model rotating discs for observed galaxies 



lAU name 


Vsys 


PA{°) 


in 


Brief comments 


(1) 


(2) 


(3) 


(4) 


(5) 


J0113+0052 


1173+ 1 


199+5 


40* 


UGC 772. S component is detached. Probable minor merger 


0122+0743W 


2953+ 2 


295+5 


37+4 


UGC 993W, component of major merger in contact 


0122+0743E 


2927+ 2 


285+5 


69 


UGC 993E, component of major merger in contact 


0335-052W 


4038+ 2 


64+8 


37+9 


centre between two knots, W comp. of major merger on HI data 


0335-052E 


4053+ 2 


53* 


37 


counter-rotation core, expanding shell, E comp. of major merger on HI data 


0822+3542 


727+ 2 


51+7 


31+15 


HI and Ha motions are decoupled 


J 1044+0353 


3865+ 1 


273+6 


51 


distorted kinematics of E-half of the disc 


1116+517 


1342+ 2 


-22+7 


50 


lack of regular rotation, probable major merging remnant 


1159+545 


3592+ 1 


-13+5 


38 


probable merger 


2236+1344N 


6176+ 2 


165+30 


31 


component of major merger in contact 


2236+1344S 


6164+ 1 


150+10 


35 


component of major merger in contact 


J012544+075957 


3012+ 3 


324+ 7 


64 


probable companion of UGC 993 


0822+3545 


742+ 4 


105* 


63* 


companion of HS 0822+3542. 



(*) — values taken from published HI maps. 

N-body simulations are performed to allow both qualitative 
and quantitative (with more detailed data) classification of 
observed velocity field s to originate due to strong interac- 
tions or mergers (e.g.. iKronberger et al.|[2006l : Ijesseit et al.l 
[20071 ). There is also substantial progress in morphology anal- 
ysis to assign ga laxies to merger products (|Conselicdl2003l : 
iLotz et aUbOOSl . and references therein). 

The main goal of these FPI Ha study was to search 
for the possible appearances of unusual, highly disturbed, 
asymmetric velocity fields in the ionized gas kinematics. The 
latter might indicate sufficiently strong interactions or var- 
ious stages of minor or major mergers. This would allow us 
to make the first step to quantify the role of interactions in 



nus model) and in the sixth panel the velocity dispersion is 
shown in colour. For both panels the intensity of Ha-line is 
superimposed by contours. The radial dependence of Vrot for 
the 'best-fitting' models of the program galaxies is presented 
in Fig. [S] 



4 DISCUSSION 

The issue of strong interactions of star-forming galax- 
ies in the aspect of their evolution i s actively studied 
in recent times (e.g. JPuech et al] I2OO6I : lYang et al.1 12OO8I : 
iDi Matteo et al.]|2008l . and references therein). Also, many 
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Figure 1. Typical individual Ha-line profiles for each spatial res- 
olution element, superimposed on the grey-scale brightness dis- 
tribution in Ho- line. The figure demonstrates the range of S/N 
ratios and of the line widths. Top panel: the region NW of 
the brightest region in galaxy SBS 0335-052E. Broadened and 
double-componen t structures of Ha- line profiles are apparent in 
the regions where llzotov et al found the evidence for out- 

flows from the similar line structures. However, our velocity res- 
olution is somewhat insufficient to derive the reliable parameters 
of the second component. Middle panel:almost the full region 
of galaxy SBS 0335-052W. The hints of the line broadening are 
seen in the outer LSB regions. Bottom panel: the NW part 
of galaxy SBS 1116-1-517. There are also indications on the line 
broadening at the outermost parts of the nebulosity. 



active SF of XMD galaxies. Of course, in some cases, the 
use of only Ha kinematics data can be insufficient to dis- 
tinguish recent interactions from other types of kinematic 
disturbance. However, the combination of ionized gas kine- 
matics with various morphology indicators, and the use of 
HI maps (when they are available) , or even of integrated HI 
profiles IConselicd |2006| ) . can lead to the correct classifica- 
tion. 

A complementary study of a subsample of XMD galax- 
ies with similar goals was conducted based on GMRT HI 
mappi ng, an d the first results a re pre s ented in Ekta et al] 
(200^, l2008l ). IChengalur eFal] ||2006D . lEkta et al.l (|2009t ). 
Some earlier results of HI mapping of XMD starb u rsting 
galaxies with VL A were presented by van Zee et al.l (|l998h 
for I Zw 18 and iPustilnik etHI (|2001bl ) for SBS 0335- 
052E,W. The Ha kinematics of ionized g as in the proto- 
t ype X MD galaxy I Zw 18 was analysed bv lPetrosian et aU 
(il995h . All these data evidence for importance of interac- 
tions in the current starbursts of studied XMD galaxies. 

If we consider the star-forming or blue compact galax- 
ies (BCGs) in general, there are very few 2D studies of Ha 
kinematics. However, for the majority of them, the morphol- 
ogy data indicate strong disturbances in the outer parts. For 
some of them, such as II Zw 40, there is clear evidence of 
a recent merger. S everal luminous BCGs were studied by 
lOstlin et al.l l|l999l ) via imaging in several bands and through 
Ha kinematics with FPL These authors concluded that all 
of their galaxies are best treated as being in the stage of 
merging. 

It is worth noting, that in such type of studies, the 
knowledge of the gas kinematics in the outer parts of a 
disc can be crucial. In particular, a good illustration is 
the luminous BCG Arp 212 (HI Zw 102). For this galaxy, 
iGarcfa-Lorenzo et al.l (|2008l ) obtained the ionized gas veloc- 
ity field for the circumnuclear region from integral-field spec- 
troscopy with INTEGRAL fiber-based system. From these 
spatial limited data, they suggested that the galaxy "shows 
a velocity field resembling a rotational system". However, 
the large-scale Ha velocity field taken with FPI, revealed a 
more complex picture where the outer HII regions belong to 
a warped pola r ring, formed via the external gas accretion 
(Moiseevl l200^V We emphasize that the spatial coverage in 
our study for all program galaxies, due to the large SCOR- 
PIO field-of-view, was sufficient to include their distant pe- 
riphery. Thus we do not miss any substantial gas velocity 
field deviations, as long as they have Ha emission above the 
detection threshold. 

First, we discuss the analysis and model fitting of Ha 
for each XMD galaxy individually (along with available HI 
maps for some of them), and then, we summarise our con- 
clusions for the subsample as a whole. 



4.1 UGC 772 = SDSS J0113+0052 

This galaxy has quite a large extent (~ 1 arcmin) and very 
complex morphology and kinematics. Observations were ob- 
tained at rather poor seeing of ~4 arcsec, but this was 
sufficient for general characterisation and modelling of its 
velocity field. The seeing probably has somewhat affected 
the picture for the 'S knot', which has a total extent of 
~ 20 X 10 arcsec. The optical light, both in broad-band con- 
tinuum and in Ha, shows clearly two distinct components. 
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The larger one (or the main body) with the overaU extent 
of ~ 50 X 20 arcsec is elongated at PA~60° and contains 
three HII regions, each having so me elongation and/or sub- 
structure. Ilzotov fc ThuanI l|2007^ obtained spectra for 3 HII 
regions in the main body on one slit. Their values of 0/H 
in all three regions are consistent to each other within the 
cited errors and can be treated as O/H for the whole main 
body. We estimated that the weighted mean for those three 
O/H values corresponds to 12-|-log(O/H)=7.29±0.05. The 
general velocity gradient in the main body is along the NS 
direction, that is almost perpendicular to the major axis of 
symmetry. 

The centre of the second distinct component (hereafter, 
UGC 772 S) is situated ~35 arcsec (or ~2.8 kpc at the 
adopted distance of 16.3 Mpc) south of the geometrical cen- 
tre of the main body. This is also elongated in both contin- 
uum and Ha, roughly in the NS direction and has a clear 
structure. There are no published measurements of O/H for 
this region. Our recent spectrum of this region at the SAO 
6-m telescope (Pustilnik et al., in preparation) results in 
12-|-log(0/H) =7.38±0.07. Despite this O/H being some- 
what larger than that in the main body, the errors of both 
values are too large to claim a real difference. The general 
velocity gradient along the EW direction (about minor axis) 
is clearly seen in this component. 

The velocity dispersion is rather low (10-20 km s~^) in 
all Hll-regions around current SF sites, where the intensity 
of Ha show local peaks. The latter is common for the ma- 
jority of other studied SF galaxies. It reaches values of 40- 
50 km s~^ in several low-brightness small periphery knots, 
where this enhancement is probably related to shocks in ex- 
panding shells. 

To zeroth order, the tilted-ring model is a reasonable 
approximation for the global motions in this object over the 
full range of the observed velocities. However, since the dis- 
tribution of Ha emission has a huge gap between the main 
body and the S component, it is difficult to find a good 
fit based only on our FP I data. Th e GMR T HI velocity 
field for this object from lEkta et al.l (|2008l ). which shows 
consistent features with our Ha velocity field in the over- 
lapping regions, is helpful for the choice of the best-fitting 
model. Therefore we accepted the value of disc inclination 
from those HI maps. The next parameters - the systemic ve- 
locity of kinematic centre, and the PA of major kinematic 
axis were derived from our Ha data. There is no well de- 
fined photometric centre. Therefore, the rotation centre was 
determined only from the symmetry of the velocity field. Its 
position also agrees (within a few arcsec) with the centre of 
HI distribution and rotation as it follows from the GMRT 
data. A model was constructed only for r > 12 arcsec; since 
in the inner region the number of points with the measured 
velocities is small and we fail to construct the stable model. 
In the outer regions the circular rotation model fits well the 
observed velocity field. The values of residual velocities in 
general do not exceed 5-10 km s~^, reaching ±15 km s^^ 
in several regions. The latter are still smaller than the ob- 
served amplitude of the rotation curve (18 km s-i). The 
largest non-random residual velocities (both negative and 
positive) are seen in the region of UGC 772 S. This looks 
like a more or less systematic gradient with the full range of 
~30 km s~\ 

Due to the patchy structure of ionized gas and the ev- 



ident problems with a simple modelization of its global ve- 
locity field, it is important to compare the Ha velocity field 
with a lower resolution, but much better sampled motions 
of HI as seen, e .g., in HI maps obtained with GMRT in 
lEkta et al.l (|2008l ). In particular, as well seen in their Fig. 9 
(right) for the FWHM beam of 15x10 arcsec, the gas veloc- 
ities observed in HI, in both the main body and UGC 772 S, 
correspond well to those measured in Ha. The general di- 
rection of the HI velocity gradient is close to NS within the 
main optical body in accord with the ionized gas motions. 
However, one can see strong deviations of HI velocities near 
UGC 772 S (and the related Hl-density peak), and in the 
southern part of the HI 'disc'. They look like an additional, 
smaller velocity component with the gradient being roughly 
along the EW direction. The other indications on the strong 
tidal HI protrusions are seen in their map in Fig. 9 (left). 
The Ha velocity field in UGC 772 S also shows a gradient 
close to the EW direction, and thus, ionized gas motions 
are coherent with those of HI. These kinematically and spa- 
tially decoupled components of UGC 772 suggest that we 
are witnessing a dwarf galaxy merger stage (the main body 
of UGC 772 and UGC 772 S), probably close to coalescence. 
An alternative interpretation of the region UGC 772 S as a 
supergiant HII region on the edge of a dwarf galaxy looks 
quite improbable due to three facts: (i) the strong distur- 
bance of the overall HI morphology, (ii) the occurrence of 
the galaxy's highest HI density peak in the edge region is 
very unusual for an isolated galaxy, and (iii) the velocity 
field in this region is clearly detached from the main galaxy 
motions. 

UG C 772 is situated in a small group of gas-rich galax- 
ies (e.g.. lEkta et al.ll2008l V Therefore, such galaxy collision 
looks quite probable. The current starburst in UGC 772 and 
in this smaller southern 'companion' galaxy is then most 
probably triggered by their collision. 

4.2 UGC 993 = HS 0122+0743 

This is one more rather extended object with very complex 
morphology and kinematics. The body, in both continuum 
and Ha, represents a group of seven bright knots (SF re- 
gions) superimposed on a LSB component. They are visually 
joined into two aggregates of similar size. The velocity field 
looks rather regular with two regions of clear gradient with 
the same direction and sign and the velocity ranges of ~100 
and ~60 km s"^ The regions roughly correspond to the W 
and E parts of the body. The velocity dispersion in most of 
the galaxy body, and in particular, around the SF knots, is 
quite low: ~10-30 km s~^. As in the previous galaxy, it is 
enhanced up to 50-60 km s^^ in several areas between HII 
regions. 

Guided by its morphology and the appearance of its ve- 
locity field, we model this object with two rotating discs in 
contact (see Figure 0). The first component, on the west- 
ern edge, is fitted with V;ys=2953±2 km s~\ PA=295°±5°, 
and with model-derived i — 37° ± 4° . Maximal rotation ve- 
locities (after inclination correction) reach ~30 km s~^. The 
second component, at the eastern edge, is best modelled with 
V;ys=2927±2 km s"\ PA=285°±3°, and i=69° (estimated 
from its axial ratio). For this component, the inclination cor- 
rected, maximal rotation velocity reaches --^50 km s~^. Re- 
sulting residual velocities appear rather small (<8 km s~^) 
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in most of the mapped area. The largest positive values of 
residuals are seen south of the bright HII region of the west- 
ern component and on the NW of the eastern component. 
They may be the tracers of current interaction or outflows. 
The two prominent, but relatively compact regions of nega- 
tive residuals (near the kinematic centre and near the SE 
edge) in the eastern component can be related to shells 
around recent starbursts. Or alternatively, the latter can 
be the tidally disturbed gas with counter-motion relative to 
that at the NW edge. For the case of modelling UGC 993 by 
a single rotating disc (as the first guess), the overall residual 
velocities were by a factor of two and more larger. 

The recent HI maps obtained at GMRT (Ekta et al. 
2009, in preparation) show clear evidence of disturbed struc- 
ture, with the large low-density HI plume stretching to the 
angular distance of ~ 1 arcmin to the south of a denser 
HI envelope that covers the optical body of the galaxy. The 
whole HI gas shows very complex kinematics. In the regions, 
where the Ha emission is well traced, HI displays the veloc- 
ity pattern similar to that presented here in Ha. 

All available morphology and kinematics data suggest 
that in UGC 993, we are witnessing a rather rare case of 
dwarf galaxy encounter in contact, when they have well 
traced regular velocity fields of both components. The kine- 
matic data suggests that the real distance between the col- 
liding galaxies is significantly larger, and their visible con- 
tact is due to projection effects. The strong SF activity in 
many HII regions, presumably triggered by the recent strong 
interaction, is spread over the components' bodies. 

From the model fit of the observed velocity field, the 
deprojected rotation velocity in the E component is a factor 
of ~1.6 larger, and hence, its estimated total mass should be 
also significantly larger. If the TuUy-Fisher scaling is valid 
for the two discs in collision (MocV"^'^), the E component 
can be ~5 times more massive. This would imply a case 
of minor merger. However, rather similar sizes of the two 
galaxies in question (as best followed namely through Ha 
velocity field) and a strong HI tidal plum mentioned above, 
are strong evidence for a major merger. The above large 
difference in rotation velocities can be due to two factors: 
(i) not properly correcting for inclination angle in the W 
component, (ii) the velocity fields themselves are affected 
by strong interactions, so the TuUy-Fisher relation may be 
not applicable. To construct a more self-consistent model, 
one needs to combine both Ha and HI velocity fields, as 
well as some additional photometry including near-IR data. 
This will be attempted in a forthcoming analysis. 

4.3 The system SBS 0335-052E,W 

The SF galaxy SBS 0335-052W is the most metal-poor gas- 
rich object with an O/ H range in different knots between 6.9 
< 12-hlog(0/H) <7.22 lllzotov et al.ll2009l). With SB S 0335- 
052E (at 22 kpc in projection. IPustilnik et aLlbOOlbl l it com- 
prises a unique pair of interacting/merging gas-rich dwarfs 
with the lowest metallicities. Both galaxies show unusually 
blue ste llar population in the outer parts of their optical 
images (jPustilnik et al.ll2004bh . The ratio M(HI)/Lb~5 for 
SBS 0335-052W is one of the greatest known. The recent 
detailed study of its disturbed H I kinematic s and m orphol- 
ogy with GMRT is presented bv lEkta etahl (120091 '). It sug- 
gests a merger of extended gas bodies soon after the first en- 



counter. However, the denser central regions of both galax- 
ies, with regions of current SF and surrounding ionized gas, 
are separated at a significant distance, and can be thus less 
susceptible to an external disturbance. Therefore, one can 
hope that kinematics of their innermost ionized gas can keep 
information on galaxy original rotation. Due to the recent 
close encounter between E and W galaxies, the gas motions 
are significantly affected by the tidal torques, especially in 
the outer parts of colliding galaxies. On the other hand, the 
gas kinematics in the innermost parts can be affected by 
current/recent SF episodes. 

One of the important issues in understanding this un- 
usual interacting pair is the following: why do the optical 
luminosities and SFRs of E and W galaxies differ by an or- 
der of magnitude, while their global HI properties are pretty 
similar? In principle, the substantial part of the answer can 
be related to differences in the geometry of collision for two 
components and the related strength of tidal action. Indeed, 
the high resolution GMRT maps of the innermost HI gas 
give evidence for a significant difference in the main veloc- 
ity gradients in E and W galaxies. For the orbital plane sky 
projection close to the EW direction, the collision of the E 
galaxy is closer to a retrograde one. Meanwhile, the inner- 
most velocity field of the W galaxy, having its main gradient 
in the N-S direction, corresponds to an almost polar colli- 
sion. This would imply a less effective tidal perturbation. 
Our FPI Ha velocity field allows a consistent check of the 
HI data. 



4.3.1 SBS 0335-05 2 W 

The observed ionized gas velocity field of the W galaxy looks 
like a combination of two components. One, centred at the 
bright HII region, shows a very small gradient, with an am- 
plitude of less than ~10 km s^^ and the direction close to 
N-S. The other refiects the gas motion between two peaks of 
Ha emission (mainly in E-W direction) and has the ampli- 
tude of ~20-25 km s~^. The 'best-fitting' model of rotating 
disc describes primarily the latter motion, which seems to 
be related to the tidal flows seen in HI maps. However, the 
PA r^lO" , accounts for the former component. Thus, the Ha 
kinematics gives additional evidence on the presence of a ve- 
locity component roughly in N-S direction, well seen in the 
high-resolution HI maps. The fainter HII region, 3 arcsec 
east of the main star-forming region, is situated along the 
direction, which is closer to the direc tion of the compac t HI 
tidal protrusion, well seen in Fig. 5 of lEkta et ahl (|2009l ). By 
analogy with the discussed opportunity of stimulated SF in 
asymmetric shells prompted by tidal outflows in the case of 
E galaxy (Ekta et al.i i20l]9l ). the similar situation can take 
place for W galaxy (albeit, with smaller power). 

4.3.2 SBS 0335-05 2E 

T he ionized gas kinem atics of SBS 0335-052E were studied 
bv llzotov et al.l (|2006l ) on the limited region of 11 x 7 arcsec 
around the central star-forming regions. In that work the 
general velocity field was not discussed. Instead, they find 
the two-component Ha-line profiles in several periphery re- 
gions, with component separation of ~50 km s~^, which 
were treated as the appearance of fast shells. In our obser- 
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vations we constructed the Ha velocity field in the region of 
~ 17 X 17 arcsec in size. 

Similar to the W component, we compare the motions 
of ionized gas in the E galaxy with that of HI gas in the 
same region. This comparison is conducted on the GMRT 
HI maps with the angular resolutions of ~20 arcsec and 9 
ar csec. Since t he res pective analysis was already performed 
m lEkta et al.l (|2009t ). this allows us to look at the observed 
Ha velocities in a more general context. Namely, there are 
two systems in the HI velocity field, which are identified 
with gas motions in elongated tidal tails. A longer and more 
rarefied one is directed in approximately SW-NE direction 
with velocities growing to the NE edge, that is in direc- 
tion opposite to the position of the companion. A shorter, 
curved denser tail starts roughly in direction to NW and 
then bends to the W companion. For this tail (protrusion), 
the velocity decreases while gas travels further from the E 
galaxy. This direction is close to that of the asymmetric 
arc of ionized gas well seen in the HST V^-band image of 
iThuan. Izotov fc Lipovetskvl l|l997r ). 

The coarse comparison of the observed Ha velocity field 
in Fig. [2] with that of the described HI one, shows qualitative 
agreement. The region of minimal line-of-sight velocities at 
~6-7 arcsec NW from the central starburst is close to the 
position of the arc. Thus, the effect of the expanding shell 
(visible as relative negative velocities on its front side) can 
also contribute to the structure of the velocity field. The 
significant swing of the isovelocity contours in the central 
region (r ^ 10 arcsec from the centre) seems may be re- 
lated to the appearance of an expanding shell. To check 
how the disturbed velocity field looks like in the presence 
of a shell (in the zeroth approximation) we consider in the 
Appendix[X] several simulated velocity fields of rotating disc 
with a shell in different positions. The observed velocity field 
of SBS 0335-052E resembles that in our Model 3, which im- 
plies an off-centre location of the shell centre on the disc 
minor axis. Based on conclusions, obtained in Appendix 1X1 
we fixed the position angle of the model disc (PA = 53°) in 
accord with the orientation of the large-scale HI structure. 
This model of rotation takes into account the main com- 
ponent of the velocity gradient (roughly from SW to NE). 
In the residual map one can clearly separate the same ve- 
locity minimum of the irregular form NW from the centre, 
related to the arc. The latter is natural to attribute to the 
approaching part of a shell, produced during a relatively re- 
cent episode of SF. Having in mind its irregular form, one 
can think of a superposition of two spatially close shells. 
Furthermore, the negative gradient in the NW direction is 
still visible, consistent with the similar fiow seen along the 
HI NW tail. 

It is worth noting the unusual behaviour in the model 
fit of rotation velocity in the very centre of the galaxy, at 
r < 2 arcsec. This can be treated as an appearance of kine- 
matically decoupled counter-rotation, which could be a trace 
of merging gaseous clouds with the specific orientations of 
their angular momenta. However, the observed rotation ve- 
locities are too small. Therefore, the more probable option 
is that this a feature of the velocity field appears due to 
non-circular motions produced by a shell formed near the 
disc centre. 

The galaxy SBS 0335-052E, unlike SBS 0335-052W, 
shows much more active star formation during the last sev- 



eral hundred Myr. This induces both large-scale shells/arcs 
and increases the amplitude of chaotic (turbulent) motions. 
The local peculiarities of Ha kinematics ar e better ac- 
quired by high-resolution spectra with VLT in llzotov et al.l 
(|2006D . In particular, they detected wide-spread regions out- 
side the central brightest clusters 1 and 2, where the line 
profiles are split on two components with velocity differ- 
ences of -50 km s-^ The 'broad' lines with FWHM up to 

— 100 km s"'^ are present in these regions as well. While 
our angular resolution is about a factor of 1.7 worse, we 
also see in the corresponding regions the enhanced veloc- 
ity dispersion of ~60-70 km s~^, corresponding to FWHM 

— 150 km s~^. Accounting for smoothing due to the worse 
seeing, these values look consistent. In fact, as the illustra- 
tion in Fig. [T] shows, our data also hint at the two-peak 
structure of line profiles in the mentioned regions, but due 
to insufficient velocity resolution they can be treated only 
as an indication. 

Summarising the overall ionized gas kinematics of the 
E galaxy, we con clude, that in accord with HI kinematics 
(|Ekta et al.ll2009l ), its direction of rotation is closer to the 
orbital plane. In general, we assume that before the merging 
galaxies experience coalescence, their innermost regions keep 
the information on their original rotation. In difference with 
the W galaxy case (polar collision), the collision of the E 
gal axy is closer to p rograde or retrograde. As emphasized 
bv lEkta et all ||2009D . the latter may be one of the reasons 
of large differences in the current and earlier SFRs between 
E and W galaxies. 

4.4 HS 0822+3542 

The Ha velocity field for this BCG is rather well fitted by 
the model of flat rotating disc. The rotation velocity curve 
raises up to —13 km s~^ at the radius r = 3 arcsec (—200 
pc) and then falls to — 7 — 9 km s~^ at radius of r — 5 — 6 
arcsec. The region with the maximal residual velocity (of 

— 10 km s~^) is positioned at a secondary peak of Ha emis- 
sion at — 7 — 8 arcsec NW of the main starburst region. This 
could be the brightest fragment of a large loop, a probable 
relic of a somewhat earlier star-forming episode. The com- 
parison of Hq kinematics with HI morphology and kinemat- 
ics (see IChengalur et al.l |2006| ) leads to interesting conclu- 
sions. At the highest angular resolution of — 6 arcsec, the 
related HI cloud is strongly elongated in the N-S direction. 
The brightest HII region is situated between two HI den- 
sity peaks. Due to low S/N ratio, the HI velocity fleld is 
constructed only for the low-resolution datacube (beam— 25 
arcsec). This has a general (but small) gradient in approxi- 
mately N-S direction. However, in the vicinity of the optical 
body one can see the swing of isovelocity line, roughly in the 
direction close to the PA of Ha velocity gradient. 

The most probable reason for the formation of decou- 
pled kinematic patterns - the internal as visible in ionized 
gas, and near polar orbits as outlined by HI - is the interac- 
tion with a nearby LSB dwarf SAO 0822-1-3545. The latter 
is also mapped in Ha in the same observations and is dis- 
cussed in a Subsection 14.91 This LSBD is situated at the 
projected distance of 3.8 arcmin (—15 kpc) and has line- 
of-sight velocity of only —15 km s^^ larger. It is —3 times 
more massive in HI a nd —0.35 mag more luminous in _B-band 
than HS 0822-1-3542 (jPustilnik et al.ll2003al : lchengalur et all 
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|2006| ). Therefore, one expects that the past interaction of the 
BCG progenitor with this LSBD could induce a starburst in 
HS 0822+3542. 

ICorbin et al.l(|2005l ). from their Hubble Space Telescope 
imaging of HS 0822+3542 with the spatial resolution of 
~5 pc, found that its central SF region consists of two 
distinct components with projected separation of ~80 pc. 
Based on this morphology, they suggest that this very com- 
pact dwarf galaxy is currently being assembled from two tiny 
unrelated components. Our Ha kinematics of this BCG, in 
difference to several other galaxies in this study, is rather 
well fit by a single disc [albeit disturbed by the earlier star- 
burst(s) and related shell(s)]. This is difficult to match with 
the idea that the BCG is currently forming from two unre- 
lated fragments. 

It is worth noting that this Ha kinematics and HI data 
from lChengalur et all l|2006l ) for HS 0822+3542 do not pro- 
vide a comprehensive picture of the overall gas motions. In- 
deed, on one hand, the regular Ha rotation within the whole 
region of optical emission (with the total extent of ~16 arc- 
sec, or ~1 kpc), has an amplitude of ~13 km s~^. On the 
other hand, HI data shows rotation with an amplitude of 
less than 5 km s~^. At the moment it is unclear how much 
the beam-smearing affects the HI velocity data. 

4.5 SDSS J1044+0353 

In the SDSS colour image, this galaxy is elongated along the 
E-W direction, with an extent of ~10 arcsec. The extent 
along the N-S direction is ~4 arcsec. The brightest part, 
the 'comet head' on the W edge, is very blue; while the 
more diffuse tail with a second, eastern 'peak' at ~4 arcsec 
is somewhat redder. The anal ysis of surface brightnes s and 
colour radial dependencies by IPapaderos et aL 1 (120081 ) con- 
firms this impression. They use the blue colours of a lower SB 
part of this object to infer the small ages of the underlying 
LSB disc. However, reservation are made on the possible sig- 
nificant contribution of nebular emission to its colour, that 
should be accounted for. Indeed, from our rather deep Hq 
image, the extent of the galaxy nebular emission (~ 12 x 8 
arcsec) appears to be similar to that of the SDSS continuum 
images. 

The two-nuclei morphology of this galaxy is sim- 
ilar to that of other studied objects (HS 2236+1344, 
SBS 1159+545), for which their Ha kinematics show good 
evidence for a merger event. The simple disc model that best 
fits the observed Ha velocity field has the following param- 
eters: VsyB = 3865+1 km s-\ PA=273°, i=51°. The bright 
rotating core is well fit by a rotation centre shifted to the 
east by ~0.5 arcsec from the photometric centre. The maxi- 
mal observed (not corrected for inclination) rotation velocity 
is 8-9 km s~^. This peak is well seen on the model rotation 
curve in Fig. [5] at r 2 arcsec. The residual velocities in 
the 'core' region (r < 2.5 arcsec) are small (< 5 km s~^). In 
outer regions the regular rotation falls significantly, and the 
residual velocities amount up to 10-15 km s~^. The fitting 
by a flat disc model with the maximal rotation velocity of 
~10 km s~^ results in rather small overall residuals. How- 
ever, the two outer regions of the E component show resid- 
uals of up to ±20 km s~^. Namely in these regions we see 
a maximal velocity dispersion of up to 50 km s~^. Besides, 
the radial dependence of rotation velocity is quite atypical 



with a substantial fall after the maximum at r ~ 2 arcsec. 
This suggests that the gas kinematics of the eastern half of 
the galaxy are signiflcantly disturbed. On the colour SDSS 
image, one can also see a redder, diffuse extension south of 
the main elongated body. All these properties are evidence 
for an unrelaxed state of this galaxy. 

The available data are not sufficient to make a confi- 
dent conclusion on the nature of this object. At the current 
level of understanding of its morphology and gas kinematics, 
the data do not contradict an idea that this is the result of 
a relatively minor, almost completed merger. Other inter- 
pretations look more problematical. The eastern component 
has a significantly smaller mass than the western one, in dif- 
ference with, e.g., the case of HS 2236+1344. This results in 
a relatively weak velocity disturbance in the W component 
and a much fainter luminosity of E component in respect of 
that for W one. High-resolution HI mapping could probably 
produce a more complete picture of its dynamics during the 
last few hundred Myr. 

4.6 SBS 1116+517 = Arp's galaxy 

This XMD BCG has a complex morphology in the central 
part, which consists of three knots with very different fiuxes 
situated in a r egion with a total extent of ~ 5 arcsec (as 
first noticed bv lArdri965l ). The outer parts also show rather 
irregular morphology with finger-like protrusions, in particu- 
lar on the S and E edges. While the full range of line-of-sight 
velocities in Ha is more than 80 km s~^, the velocity field 
is rather complex. It is hard to see tracers of regular rota- 
tion. Nevertheless, we try to construct the model of rotating 
disc for radial distances r < 9 arcsec, which we extrapolate 
to the outer region. The centre of the inner isophotes was 
set as the centre of rotation. As one expects, the rotation 
is indeed small: its maximal observed line-of-sight veloci- 
ties are of order 5-10 km s~^, while the residual velocities 
are ~20 km s~^ in the central part and reach the value of 
~50 km s~^ at the E and SE periphery. The region with the 
most peculiar (redshifted) velocities is situated South-East 
from the bright 'central' knot and possibly represents the 
result of a merger (or infall to the disc). 

While the 'best-fitting' model with rotation can account 
for maximal rotation velocity of ~10-12 km s^^, the result- 
ing residual velocity field shows the regions with values down 
to -20 km s~^ and up to +50 km s~^. This very disturbed 
velocity field, especially in outer parts of an ionized gas body, 
suggests strong perturbation. The presence of expanding 
shells from recent starbursts is probably one of the possible 
reasons. If one assigns the regions with the most negative 
velocities on the residual map to approaching fronts of such 
shells, their expansion velocities are ~20 km s~^ and the 
total sizes are ~3-5 arcsec (or ~300-500 pc at the adopted 
distance of 23 Mpc). The sizes are quite typical of ionized-gas 
shells, visible in similar galaxies. However, the large positive 
residual velocity at the eastern edge of the galaxy is diffi- 
cult to reconcile with shells. The morphology of Ha-emission 
region and velocity fields, both the observed and residual, 
suggests strong external disturbance. No candidate galaxy 
capable of producing that strong of a disturbance is visible 
in the BCG environment. Therefore, the most natural in- 
terpretation that emerges from the study of SBS 1116+517 
morphology and gas kinematics is the recent major merger 
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of dwarf galaxies, since for a minor merger one expects a rel- 
atively weak effect on the overall velocity field of the more 
massive component. If this picture is correct, one should ex- 
pect that the HI morphology and velocity field are highly 
disturbed. Thus, this galaxy requires follow-up HI mapping 
with angular resolution of 3-5 arcsec. 



4.7 SBS 1159+545 

On the SDSS image SBS 1159-1-545 consists of two very blue 
knots (separated by ~ 5 arcsec, or ~1.3 kpc at the adopted 
distance of 52.5 Mpc), elongated roughly along NE-SW di- 
rection, with the luminosity difference of a factor ~2. There 
are also redder lower brightness fragments; one is situated 
between the two blue knots and the other, with the total ex- 
tent of ~ 5 arcsec, stretches as a tail of the fainter, blue SW 
knot. The SW knot itself has an extended structure, roughly 
oriented along E-W, and its redder tail looks like a contin- 
uation of this blue knot. On the other hand, the brighter 
NE knot also looks elongated on our narrow red continuum 
image, in a direction of S-N. This overall complex, curved 
morphology indicates a non-equilibrium state of this object. 

The full range of Ha velocities seen in our data is 
~40 km s~^. The velocity field is complex, and does not re- 
semble a projected rotation. Moreover, it resembles that of 
some other objects in this study, with 'two-nuclei' morphol- 
ogy, which we have argued are mergers close to the full co- 
alescence (SDSS J1044-H0353 and HS 2236-H344). Namely, 
there are two local maxima in the line-of-sight velocity, close 
to the positions of the Ha-emission peaks and a clear min- 
imum along the border, or 'middle' lane, between these SF 
knots. From the position of this middle lane, the observed 
velocities grow to both NE and SW directions. However, af- 
ter reaching the local maxima near the Ha emission peaks 
they again fall. The largest velocities are seen on the W and 
E edges of the NE knot as well as at similar outer parts of 
the SW knot. 

The velocity field model was built for the range r < 5 
arcsec (with the centre at the NE knot) and was extrapo- 
lated to the SW knot. The 'best-fitting' model for the fiat 
rotating disc centred on the brighter NE knot can account 
for only an amplitude of ~8 km s~^ along the line-of-sight. 
The residual velocities are significantly larger, amounting 
on the knot's east edge of 15-20 km s"'^. Similar indica- 
tions of an overdisturbed velocity field come from the map 
of velocity dispersion. While in the regions around bright 
knots, this parameter does not exceed 15-20 km s~^. In the 
E part of NE knot and the W part of the SW knot, this 
exceeds 30-50 km s"^ The PA=-13° of the 'best-fitting' 
kinematic model significantly differs from PAk^ 30° of the 
apparent 'disc' plane - if the two knots were HII regions in 
the same disc galaxy. But this is consistent with the orien- 
tation of the NE knot itself, as seen from our narrow band 
continuum image. An attempt to interpret the whole veloc- 
ity field as a result of combination of two expanding shells 
does not work. In this model, it is difficult to understand 
how the gas in the middle lane acquired its low observed 
velocity, while the line-of-sight velocities in the directions 
of expected centres of shells (NE and SW knots) show local 
maxima. Summarising this velocity field analysis, we suggest 
that from both the BCG morphology and gas velocity field. 



the most plausible interpretation of the data is the merger 
of two dwarf galaxies. 



4.8 HS 2236+1344 

This XMD BCG also has the two-nuclei morphology with 
disturbed outer parts (Fig. |4]). The separation of two bright 
knots, oriented in the N-S direction, is ~ 3 arcsec (or ~1.2 
kpc in projection at the adopted distance of 85 Mpc). As 
seen in deep images (not shown here), the latter includes 
small tails protruding from both the northern edge of the N 
component (to the W) and from the southern part of the S 
component (Pustilnik et al., in preparation). These features 
already hint on the tid al perturbations in both components. 
iPramskii et al.l l|2003l ) studied the Ha velocity field of this 
object with a long slit, positioned along the N-S direction. 
They discovered in the position-velocity diagram that the 
line-of-sight velocity makes a jump in the middle between 
the two knots. Our Ha velocity field confirms that finding 
and adds important new features, including high velocity 
dispersion and large residual velocities near the regions of 
current starbursts for any model of single rotating disc. 

The total velocity range is ~45 km s~^. We attempted 
to fit the whole observed velocity field with the model of 
single fiat rotating disc. Two variants of the circular motion 
model were constructed: (i) rotation centre coincides with 
the brighter (Southern) knot; (ii) rotation centre is fixed in 
the middle between Northern and Southern knots. However, 
both variants fail to reproduce the observed velocity gradi- 
ents. For an alternative approach, the approximation of the 
overall Ha velocity field by two rotating discs, positioned 
near the centres of S and N knots, gives a good fit with 
maximal rotation velocities of ~25 km s~^ in both compo- 
nents (see Tableland Fig. [3] and [S}. These discs have close 
PA, Vsys, i and maximal Vrot. 

Due to lack of space, we focus on the model with two 
discs which provides a better fit to the observations than ei- 
ther single disc model. The map of velocity dispersion shows 
this parameter to be rather large in this BCG. Even its min- 
imal values, located in the regions of both bright knots, are 
of ~25 km s^^. In the adjacent regions, this parameter raises 
till 35-40 km s^^, while in the outer parts, the velocity dis- 
persion exceeds 50 km s~^. The observed kinematics can not 
be explained as well by giant ionized shells in a single rotat- 
ing disc. The latter are commonly expected to show negative 
relative velocities in their front edges, that is near the cen- 
tres of star-forming regions, while the observed residuals are 
positive. 

Taking into account all available morphological and 
kinematics data on this BCG, we conclude that the most 
likely nature of this object is a close strongly interacting 
pair of dwarf galaxies with comparable masses. It is quite 
probable that this pair is in the stage of merger close to full 
coalescence. This also provides a natural interpretation of 
th e starbursts in both components. Indeed, as summarised 
bv lBegum et al.l l|2008l ). the late- type quiet dwarfs with that 
low Vrot have typical Mb in the range of -12.5 to -14.9. 
Therefore, the Mb=-16.9 for HS 2236+1344 indicates a very 
strong starburst. 
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4.9 Properties of two by-product LSB galaxies 

In this section we briefly discuss two galaxies that appeared 
sufficiently close on spatial and velocity separation to the 
main studied XMD galaxy and were mapped along with the 
main target. 



4.9.1 Anon J012544+075957 

In the FPI datacube for UGC 993, an LSB dwarf galaxy 
Anon J012544+075957 appeared (see bottom of Table [g]) at 
2.5 arcmin East and ~ 40 arcsec North from the main target 
(~29 kpc in projection). Its line-of-sight velocity (as found 
here from its Ha velocity field) is ~70 km s"'^ larger than the 
mean velocity of the UGC 993 system. Its Hq image (Fig. |4]) 
displays three main regions elongated roughly from NE to 
SW, forming the main body with the total extent of ~ 25 
arcsec and a separate faint knot SW from this. The veloc- 
ity range along the body is ~70 km s~^. The 'best-fitting' 
model parameters for a rotating thin disc are as follows: 
Vsys=3012 km s"\ PA=324°, i=64°. The resufiing rotation 
curve reaches a maximum value of ~25 km s~^ (Fig. [SJ at 
the edge of the Ha map {r ~ 12 arcsec). The gas morphology 
of this LSB galaxy is clearly disturbed and elongated nearly 
perpendicular to the plane of rotation. It is interesting that 
in the blue continuum, the galaxy looks less elongated than 
in Hq. However, there are two extensions approximately to 
the N and S, close to the elongation of Ha emission. The 
velocity dispersion near the peaks of Ha emission is a mod- 
est ~25 km s~^. In some of the outer regions, it raises up to 
60-70 km s^^, but this is probably attributed to a low S/N 
ratio. 

Along with the already discussed two kinematically de- 
coupled components of UGC 993, this dwarf Irr galaxy likely 
comprises a progenitor triplet of low-mass galaxies, whose 
common dynamical evolution resulted in a recent merger, 
as seen in UGC 993. The general problem of such dwarf 
bound systems was addressed, in particular, in ITuIIv et al.l 
l|2006l ). 

4.9.2 SAO 0822+3545 

Another by-product object, SAO 0822-1-3545, the LSB dwarf 
galaxy, kno wn as a companion of HS 0822-1-3542, was already 
discussed in^ Pustilnik et all (|2003al ). Its unusually blue BVR 
colours imply a rather small age for its main stellar popu- 
lation. Low level SF in this LSBD was noted in the latter 
study based on the detection of a low EW Hq emission in 
the long-slit spectrum taken along the major axis. The new 
data show the 2D distribution of Hq emission (see Fig. |4]) 
and allow us to perform a more careful analysis of SF issues. 

In particular, our Hq image shows that the current SF 
in this LSB dwarf is highly inhomogeneous. It is concen- 
trated in two close HII regions No. 1 and 2 (~ 4 arcsec, or 
~0.26 kpc in between) at the Northern edge of the main 
galaxy body, where no light concentration is seen in the 
broad-band continuum images. A fainter Ha protrusion is 
seen in the W part of the faint irregular structure. The Hq 
kinematics are complex. On one hand, due to very limited 
spatial coverage of Hq emission over the LSBD body, it is 
difficult to derive the parameters of a rotating disc model 



based only on Hq data. On the other hand, the Hq veloci- 
ties are in good agreement in the overlapping regions wit h 
those observed in HI as presented in lChengalur et al] (|2006l ). 
Based on this fact, we fixed the position of the rotation cen- 
tre, PA, and i, as they emerge from the mentioned HI maps. 

As usual, the velocity dispersion near the brightest Ha 
knots is rather low: ~5-15 km s~^. In several regions outside 
these two knots, where Hq is still seen - at ~ 4 arcsec SE of 
knot No. 1 and in the W protrusion - the velocity dispersion 
is large, up to 60-70 km s~^ . This could be evidence of recent 
strong disturbances of the galaxy's ISM. We also notice that 
the total Hq flux in this LSB galaxy is ~25 times smaller 
than that in its companion BCG HS 0822-f -3542, while its 
total b lue luminosity is ~30 per cent higher jPustilnik et al.l 
l2003al ). This seemingly reflects the difference between the 
SF properties in the two interacting galaxies, related to the 
difference in their density distribution and the strength of 
the tidal trigger. 



5 SUMMARY 

Summarising the results and discussion above, we draw the 
following conclusions: 

(i) The ionized gas kinematics in very low-metallicity 
star-forming galaxies, studied with FPI Ha observations are 
signiflcantly disturbed and rarely can be well fltted by only 
single disc with regular rotation. 

(ii) Several of our galaxies show more or less clear ev- 
idence from morphology and Ha velocity fields for var- 
ious stages of mergers (UGC 993, SDSS J1044-h0353, 
SBS 1116-H517, SBS 1159+545, HS 2236-H344). In two cases 
we have good evidence for two independently rotating discs 
(UGC 993 and HS 2236-f 1344). 

(iii) In the other galaxies, the rotation component of the 
overall velocity field is important, but large disturbances ap- 
pear. The residuals of the 'best-fitting' rotation model imply 
either a recent merger, or sufficiently strong disturbance by 
nearby galaxies (HS 0822-^3542 and UGC 772). 

(iv) Probable starburst-induced shells were identified in 
several galaxies (UGC 772, UGC 993, SBS 0335-052E, 
SBS 1116-1-517) through their ionized gas velocity fields and 
velocity dispersion maps. We presented the results of simple 
simulations of the expected velocity patterns which one can 
observe in dwarf galaxies with expanding shells. To first or- 
der, this allowed a feel for the effect of shells on the results 
of the tilted-ring model. 

(v) The interacting/merging nature of the binary system 
of the weU separated XMD galaxies SBS 0335-052E and W 
is best evident from their HI morphology and kinematics 
data. Despite a relatively large mutual distance, the tidal 
action of each component to the other clearly affects the 
gas dynamics of these very gas-rich objects and triggers 
the current SF burst and ve ry likely the pre vious major 
SF episode (as emphasized bv lEkta et al.l [20091 ). To under- 
stand this unique interacting system (a nice representative 
of high-redshift young galaxies) in more detail, one needs 
a wide grid of models of interacting gas- rich galaxies, like, 
e.g., "Identikit" (iBarnes fc Hibbardil2009t ) . but including SF 
processes. 

(vi) The Hq images and velocity fields for two LSB 
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dwarfs, Anon J012544+075957 and SAO 0822+3545, com- 
panions of two program XMD galaxies, are obtained and 
analysed. They can be used in statistical studies of SF in 
LSB dwarfs. The star formation in SAO 0822+3545 is highly 
asymmetric and takes place mainly in two knots at the 
northern edge, that is likely induced by the recent inter- 
action with the nearby XMD BCG HS 0822+3542. 

The statistics of XMD starbursting galaxies, for which 
kinematics of gas were studied in detail, is still insuffi- 
cient. Nevertheless, the results of our FPI study of the ion- 
ized gas kinematics in the subsample of nine XMD star- 
forming galaxies, along with the complementary results of 
the GMRT HI study of a pa r t of these and other XMD galax- 
ies l|Ekta et al.ll2008l . l2009l : lEkta fc Chengaluij|2010i ). indi- 
cate that strong interactions and mergers of very metal-poor 
dwarf galaxies are one of the major or significant factors trig- 
gering their current and recent starbursts. In particular, this 
is valid for the six most metal-poor (12+log(O/H)<7.30) 
dwarf starbursting galaxies. Merger-induced starbursts are 
consistent with the idea that the progenitors of such rare 
objects either are old and have been evolving very slowly on 
the cosmological timescale before the current starburst have 
occurred, or they are extremely metal-poor because they are 
comparatively young, and thus began their star formation 
and chemical enrichment with large delay. 
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APPENDIX A: ON THE VELOCITY FIELD OF 
DISC WITH SUPERIMPOSED SHELL 

The tilted-ring fit is very po pular for analysis of gas velocity 
fields in galactic discs (see ITeubenl ()2002i) for detailed re- 
view). In this approach, it is easy to connect the radial trends 
of kinematically determined PA and i with real changes of 
disc geometrical parameters (warp, polar rings, etc.) and/or 
with regular non-circular motions caused by spiral arms or 
bar. However, in dwarf galaxies the SF-induced expanding 
shells can significantly disturb the disc velocity field since 
the amplitude of the rotation curve is rather small and the 
size of a shell can be comparable with the size of the disc. If 
we allow for variations of kinematical PA with radius, it can 
lead to a mistaken conclusion on disc orientation or about 
the character of non-circular motions. In this Appendix we 
examine the effect of expanding shells on the tilted-ring fit- 
ting of artificial velocity fields. 

The 'observed' velocities in our two-component model 
are simply the sum of line-of-sight velocities for the disc and 
the shell: 

Vobs = Vdisc + Vshell = Vsys + Vrot COS Iy3 siu 1 + Vshell ■ ( Al) 

This approximation assumes the equal contribution of both 
components in the total velocity pattern in the overlapping 
regions, i.e., the same brightness of the components in each 
point. Of course, this is not similar to the real situation. 
However, we believe that eq. ([Aip can be used as a first 
approximation. One of the constraints used in the models 
below is that VshM should not exceed the typical FWHM 
of the observed emission lines (about 30 — 50 km s~^). Oth- 
erwise we would detect the two-component line profiles. 

The azimuthal angle in the plane of galaxy tp is con- 
nected with the position angle in the sky by the relation: 
tanip = tan(P^ — PAo)/ cos i. The rising rotation curve is 
constructed by the relation: 



_ 2 R 
Vrot — — Vmax arctaii — , 
TT a 



(A2) 



with h = 1 arcsec and the maximal velocity Vmax = 
30 km s~^ which is typical of dwarf galaxies in our sam- 
ple. 

For the shell, we use a model of a thin expanding 
hemisphere. This implies that we see only its approaching 
(blueshifted) side, whereas the redshifted one is obscured by 
an interior extinction. The model for the shell is: 



Vshell = — Vea 



{X - Xs)'^ + {y- VsY 



(A3) 



Here Xs and j/s are sky-plane coordinates of the shell centre 
relative to the disc centre. The shell has a radius of Rs = 3.5 
arcsec and the expansion velocity of Vexp =15 km s~^. 

Using eq. (|A1|) . we created models for moderately in- 
clined disc {i = 45°, PAo — 45°) with the maximal radius 
of 10 arcsec and with various values of Xs and i/s- A systemic 
velocity of Vsys = 1000 km s~^ was adopted, that is repre- 
sentative of our sample. Its value has no effect on model 
isovelocity patterns. 

The models were calculated on the grid with cell size 0.1 
arcsec. Then the constructed models were convolved with 
the 2D Gaussian with FW HM — 1.5 arcsec that simulated 
the typical seeing effect, and were rebinned with the pixel 



size of 0.5 arcsec in accord with our observational conditions. 
Fig. lAll shows examples of typical simulated velocity fields. 

Model represents the disc without shell. In Model 1 
the shell is located in the disc centre. In Model 2 and 
Model 3, the centre of the shell is shifted along the ma- 
jor and the minor axes, respectively. In Model 4 the shell 
centre is shifted in the direction of PA — 180°. 

The simulated velocity fields were analysed with the 
same tilted-ring model method under two approximations. 
Method I allows for radial variations of PA and Vsys, while 
in Method II these values were fixed and only the rotation 
velocity was a free parameter. The inclination angle and the 
rotation centre were fixed in both Methods. Fig. lAll shows 
the residual velocities in the both approximations. The re- 
sults of fitting are shown in Fig. IA2I These figures demon- 
strate how an expanding shell affects the observed velocity 
field and changes the estimates of best-fitting parameters. 
Namely, the central location of the shell, or its offset along 
the disc major axis provoke the radial variations of Vsys (see 
Models 1 and 2). The shell location beyond the disc major 
axis also gives the systematic errors in the PA estimations 
(Models 3 and 4). When the shell centre lies on the minor 
axis, the radial trend of PA reach the maximal amplitude 
(~20° in our simulations). That large variation of PA can 
result in erroneous conclusions concerning the disc structure, 
like the presence of the inner warp or even of a nuclear disc. 

The residual velocities in the case of Method I seem 
spread over a larger portion of the disc than the real area 
of the shell (Fig. IA1|) . The residuals are relatively small, 
usually about ±5-7 km s~^. Moreover, the peak of negative 
('blue') residuals usually does not coincide with the centre 
of the shell. In contrast, if PA and Vsys are constant (i.e., 
for Method II), the residual velocities pattern is in a good 
agreement with the location of the shell, and the value of 
the negative residuals is consistent with Vexp, with except 
of the special case of Model 2. 

It is worth noting that the rotation curve derived from 
the simulated velocity field can be distorted independently 
of the analysis method, as can be seen for Vrot plots for 
Models 2 and 4 in Fig. IA2I However, in the Method II ap- 
proach, we can unambiguously recognise from the residual 
map the perturbed portion of the disc and then derive the 
real rotation curve on the second iteration after masking the 
shell region. The latter is impossible if PA and Vsys are free 
parameters of the fitting, because in this case it is difficult 
to recognise the shell location. 

Based on the results presented above, we have analysed 
the velocity fields of dwarf galaxies using the 'tilted-ring' 
fitting assuming that PA and Vsys are constant with radius. 
When the first approximation was done, we masked the re- 
gions with large residuals (usually, larger than 5 — 10 km s~^) 
and repeated the fitting. The final model velocity field and 
the rotation curve, i.e., parameters of the disc orientation 
were calculated after several such iterations. 
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UGC772 UGC993 SBS0335-052W SBS0335-052E 
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Figure 2. Images and FPI maps for four program galaxies UGC 772, UGC 993, SBS 0335-052W and SBS 0335-052E. From top to 
bottom: broad-band and Ha images in logaritlimic intensity scale; the line-of-sight velocity field in Ha and model of circular-rotated 
disc, the scale is in km s"'^ relative to the systemic velocity which is corresponded to thick contour; the maps of residual velocities and of 
velocity dispersion. In two latter maps the isophotes of Ha-line intensity are also superimposed. The cross marks the kinematic centre. 
For UGC 993 and HS 2236-1-1344, centres of both discussed components are marked. The step between isovelocities is 5 km s~^, with 
the exception of UGC 993, SBS 0335-052E and Anon J012544+075957 where the step of 10 km s"! was used. The scale bar in the top 
panel corresponds to 1 kpc. 
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Figure 4. Same as in Fig. [2 but for XMD galaxy HS 2236+1344 and for two by-product LSB galaxies Anon J012544+075957 and 
SAO 0822+3545. 
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Figure 5. Radial dependence of rotation velocity in the best-fitting models for all discussed objects. 
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Figure Al. The simulated maps, (a) — the relative position of the disc (grey) and shell (dark), (b) - the simulated velocity field (after 
the subtraction of Vsys = 1000 km s~^). (c) - the residual velocity field after the subtraction of the tilted-ring model (Method I) which 
allows variations of PA and Vsys- (d) - the same residual field, but for Method II {PA, Vsys =const) 
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Figure A2. The fitting parameters for the simulated velocity fields: rotation velocity (top), position angle (middle) and systemic velocity 
(bottom). Blue points and red diamonds show results for Method I and Method II, respectively. 
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